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Abstract 

The new exact formulas for the attractive Casimir force acting on 
each of the two identical perfectly conducting plates moving freely 
inside an infinite perfectly conducting cylinder with the same cross 
section are derived at zero and finite temperatures by making use of 
the zeta function technique. The long and short distance behaviour 
of the plates' free energy is investigated. 

1 Introduction 

Recently a new geometry in the Casimir effect pQ, a piston geometry, has been 
introduced in a 2D Dirichlet model [2]. Generally the piston is located in a 
semi-infinite cylinder closed at its head. The piston is perpendicular to the 
walls of the cylinder and can move freely inside it. The cross sections of the 
piston and cylinder coincide. Physically this means that the approximation 
is valid when the distance between the piston and the walls of a cylinder is 
small in comparison with the piston size. 

In paper [3j a perfectly conducting square piston at zero temperature 
was investigated in 3D model in the electromagnetic and scalar case. The 
exact formula (Eq.(6) in [3]) for the force on a square piston was written 
in the electromagnetic case. Also the limit of short distances was found for 
arbitrary cross sections at zero temperature (Eq.(7) in [3]). 
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A dilute circular piston and cylinder were studied perturbatively in [3]. 
In this case the force on two plates inside a cylinder and the force in a piston 
geometry differ essentially. The force in a piston geometry can even change 
sign in this approximation for thin enough walls of the material. Other 
examples of pistons in a scalar case were presented in El [7] . 

In this paper we consider a 3D model in the electromagnetic case with 
physically realistic perfectly conducting boundary conditions. Two identical 
plates move freely inside an infinite cylinder with the same cross section, 
which is arbitrary. The plates are perpendicular to the walls of the cylinder. 
In Section 2 we derive the new exact zero temperature result (j2"6j) . (1281) for 
the Casimir energy of two identical parallel plates inside a cylinder with 
an arbitrary cross section by making use of the zeta function technique [8, 
[H]. Two special cases of (|2"£1) . the exact results for rectangular and circular 
cylinders, are briefly discussed. Also we discuss the important property of 
the result f l28l - its exponential damping at long distances ( 1311) . In Section 3 
the new exact result (134")) for the free energy of two identical parallel plates 
inside an infinite cylinder is derived. In the long distance limit the new high 
temperature result (I36p for the free energy is obtained. In Section 4 the heat 
kernel technique jTOJ El EG] is applied to derive the leading short distance 
behaviour of the free energy. Also we prove that in the short distance limit 
of the high temperature result (1361) the known high temperature result for 
two infinite perfectly conducting parallel plates floTj) is obtained. 

We take h = c = 1. 

2 Zero temperature results 

Our aim is to calculate the Casimir energy of interaction and the force be- 
tween the two identical parallel plates of an arbitrary cross section inside an 
infinite cylinder of the same cross section (the plates are perpendicular to 
the walls of the cylinder). 

TE and TM eigenfrequencies of the perfectly conducting cylindrical res- 
onator with an arbitrary cross section M can be written as follows. For 
TE modes (E z = 0) inside the perfectly conducting cylindrical resonator 
[0, a] x M with n being a normal to the boundary dM of an arbitrary cross 
section M the magnetic field B z (x,y, z) and eigenfrequencies ujte are deter- 
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mined by: 



+00 ^ 

B z (x,y,z)= B u sin (—)9i(x,y), (1) 

i=l,l=l 

A^g l (x,y) = -\ 2 N g l (x,y) (2) 

d9 ^ y) =0 (3) 
an om 

^Ti? = ( — ) + \jv, ,A iA r^0, Z,i = l.. + 00. (4) 

The other components of the magnetic and electric fields can be expressed 
via B z (x,y,z). 

For the TM modes (B z = 0) inside the same perfectly conducting cylin- 
drical resonator the electric field E z (x,y,z) and eigenfrequencies ujtm are 
determined by: 

+00 ^ 

E z (x,y,z)= Y E kicos( }fk(x,y), (5) 



V a 

l=0,k=l 



A^f k (x,y) = -\t D f k (x,y) (6) 
fk(x,y)\ dM = (7) 

w ™ = (^) +A ^ D ' Z = - + °° ,k = 1 - + 00 ( 8 ) 

In zeta function regularization scheme the Casimir energy is defined as 
follows: 

E = 1{^te+J2 uj tm)\ s=0 , (9) 

where the sum has to be calculated for large positive values of s and after 
that an analytical continuation to the value s = is performed. 

Alternatively one can define the Casimir energy via a zero temperature 
one loop effective action W (7\ is a time interval here): 

W = —ETi (10) 
E = -C(0) (11) 

£ £T£->H;)V+')) (12) 

U1 TE ,U) TM 
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After integration over p in ffT2l one can see that definitions @ and ffTTj) 
coincide. 

In every Casimir sum it is convenient to write: 



+00 ^ , 

^exp(-t/ 2 ) = -e 3 (0,- 



£\ 1 
2" 



(13) 



For the first term on the right-hand side of ( fT3l) we use the property of the 
theta function # 3 (0,x): 



3 (O,x) = -^(o,- 
and the value of the integral 



(14) 



dir- ! rxp( -pt- |) = 2(iy K a (2y/pq) 



(15) 



for nonzero values of n to rewrite the Neumann zeta function (n{s) (arising 
from TE modes) in the form: 



Cn(s) 



+00 



7rr((a-l)/2) /av^Tp 



2T(s/2) 

-r^72)^ a> /AS7+? 

E 



2 X 1 



7T 



+ 



\- 20T / 7T 2 Z \ , 



J Rr 2= i(2aZJA? JV +p 2 



Vir((s-l)/2) ( a\ 



N 



4ar(s/2) V 7T 



(16) 



The Neumann part of the Casimir energy is given by: 



2vr2 
+ 



a v ^ r 00 dp / 
2^7^ 2^1 

AiiV 



4v + ^ 



5=0 A*, 



(17) 



s=0 



Here we used K_ 1 / 2 (x) = \J tx j (2 x) exp(— x). 
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The Dirichlet part of the Casimir energy (from TM modes) is obtained 
by analogy: 



E D = J2 |^ln(l - exp(-2a v /A^ +p»)) + 



5>i 



— s 



;is) 



s=0 



The electromagnetic Casimir energy of a perfectly conducting resonator 
of the length a and an arbitrary cross section is given by the sum of ffTTl) and 
(EED : 



E =E / J 5^ K 1 " ex P (-2a v /A? JV +^)) + 
+ E r°°^^n(l-exp(-2a v /ALT^)) + 



s=0 



+ 



s=0 



(19) 
(20) 
(21) 
(22) 
(23) 



s=0 



The terms 



E, 



cylinder 



P 



\VL d ^ N+ 



+ 



s=0 



s=0 



(24) 



yield the electromagnetic Casimir energy for a unit length of a perfectly 
conducting infinite cylinder with the same cross section M as the resonator 
under consideration. 

For rectangular boxes it was generally believed [331 H3J that the repulsive 
contribution to the force acting on two parallel opposite sides of a single box 
(and resulting here from (12111221) ) could be measured in experiment. However, 
the terms (I21ti22|) are closely related to the Casimir energy of an infinite in z 
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Figure 1: Two plates inside an infinite cylinder 



direction cylinder when there are no plates inside (see eq. (124jl ). Imagine that 
the box is large in z direction. Its Casimir energy and the Casimir energy 
of an infinite in z direction cylinder coincide when the two opposite z sides 
of the box are located at spatial infinity. To calculate the energy change 
between these two configurations and the force on a z side of the box one 
should subtract the energy of an infinite cylinder from the expression for E 
( |T9ll23l) when the box sides are infinitely far from each other. Then the force 
on a z side of the box is equal to zero for infinite distance a between box z 
sides. 

For the experimental check of the Casimir energy one should measure 
the force somehow. One can insert two parallel perfectly conducting plates 
inside an infinite perfectly conducting cylinder and measure the force acting 
on one of the plates as it is being moved through the cylinder. The distance 
between the inserted plates is a. 

To calculate the force on each of the two plates inside a cylinder with the 
cross section M one can perform the following gedanken experiment that was 
frequently used to calculate the Casimir force between two infinite parallel 
perfectly conducting plates. Imagine that 4 parallel plates are inserted inside 
an infinite cylinder and then 2 exterior plates are moved to spatial infinity. 
This situation is exactly equivalent to 3 perfectly conducting cavities touching 
each other. From the energy of this system one has to subtract the Casimir 
energy of an infinite cylinder without plates inside it, only then do we obtain 
the energy of interaction between the interior parallel plates, the one that 
can be measured in the experiment. Doing so we obtain the attractive force 
on each interior plate inside the cylinder: 



F(a) 



d£(a) 



(25) 



da 




(26) 
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the sum here is over all TE and TM eigenfrequencies tu wave for a cylinder 
with the cross section M and an infinite length. 

In fact, the final results for this geometry should coincide with the results 
for the three plates' piston geometry when one of the three piston plates (the 
exterior plate) is moved to infinity. Also one can immediately obtain the 
result for three plates' system inside a cylinder, which is exactly the piston 
geometry, employing the same arguments and renormalization scheme. In 
the three plates' system the force on the interior plate is equal to the sum of 
the forces acting on this plate from the two exterior plates, i.e. the piston 
geometry can be effectively considered as 2 two plates' systems. It should be 
emphasized that these assertions are valid for perfect boundary conditions. 

By making use of an identity 



1 f +c °dp, , . r- - N X^KA2lXa) 

2 4r ln(l ~ ex P(~WA 2 +P 2 )) ="2z£ "V^ (27) 
J -°° i=i 

one can rewrite (126]) in the form: 

c , \ 1 V^/v^ ^kDKi{2l\ kD a) sr-^ X iN K 1 (2lX iN a)\ 

£ w = -^2^{2^ ] + \ ) ■ ( 28 ) 



The results ( 1261) , (1281) are our main zero temperature results. Our results 
are exact for an arbitrary curved geometry of a cylinder. This fact may be 
important for the experimental check of the Casimir effect in piston related 
geometries with curved boundaries. One can choose an arbitrary curved 
plate geometry M, for this geometry the eigenvalues of the two dimensional 
Dirichlet and Neumann problems XkD, Xny can be found numerically. After 
that the exact expressions ( 1261) . ( 1281) can be used to obtain the Casimir force 
on a plate. In fact, similar in the form mathematical results can be obtained 
in the case of a one dimensional massive scalar field theory on a manifold 
with boundary, however, in our case an infinite number of "masses" XkD, 
appear in the theory due to existence of the cylinder cross section. 

For convenience of the reader we write explicitly two special cases: rect- 
angular and circular cylinders. For a rectangular cylinder with the sides b 
and c the exact Casimir energy of two plates inside it can be written as: 

Srect{a) = -Y J £ \An 2 /frj + n 2 /c 2 /&2 + n2/c2 y (29) 

1=1 m,n=—oo 

The prime means that the term m = n = is omitted in the sum. 
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For a circular cylinder the eigenvalues of the two dimensional Laplace 
operator \ kD , X iN are determined by the roots of Bessel functions and deriva- 
tives of Bessel functions. The exact Casimir energy of two circular plates of 
the radius b separated by a distance a inside an infinite circular cylinder of 
the radius b is given by: 



+00 +00 

1=1 u=0 j 



1 fiDvjK^lfiDvja/b) + U NvjK^lnxyja/b) 

(30) 



Jv{V>Dvj) — 0, J u (fiNuj) — 0. 

The sum is over positive HDvj and fiNuj- 

The leading asymptotic behaviour of £{a) for long distances Ai^a ^> 
1, Auva ^> 1 is determined by the lowest positive eigenvalues of the two 
dimensional Dirichlet and Neumann problems Aid,Aiat: 

£(a)U 1D a»i,A 1Jva »i ~ - I j=(^/^e- 2X ^ a + ^^e- 2 ^, (31) 



so the Casimir force between the two plates in a cylinder is exponentially 
small for long distances. This important property of the solution is due to 
the gap in the frequency spectrum or, in other words, it is due to the finite 
size of the cross section of the cylinder. 

3 Finite temperature results 

To get the free energy J 7 (a, f$) for bosons at nonzero temperatures (T = 1/(3) 
one has to make the substitutions: 

2nm 

V ^Pm = (32) 

+00 j 1 +°° 



L 



Thus the free energy describing the interaction of the two parallel perfectly 
conducting plates inside an infinite perfectly conducting cylinder of an arbi- 
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trary cross section has the form: 



, +00 

*E E 2 ln ( 1 - exp (- 2a v A ^ + ^))+ 



^kD rn—-^ 

+ *E E 2 ln ( 1 - ex p(- 2a v A ^+^))- (34) 

This is our central finite temperature result. Note that Any 7^ 0. For rect- 
angular and circular cylinders one can substitute the explicitly known AkD, 
A iN to flSD in analogy to (|29l) and (|30l). 

The attractive force between the plates inside an infinite cylinder of the 
same cross section at nonzero temperatures is given by: 



dF{a,p) 1^ "td 

l fin /lZ^p Yn f%M m , 



9a /3 ^-^ exp(2acj TD ) — 1 



/3^exp(2aa; TiV )-l' 



\ 2 



Here u; T d = \/Pm + A^, D and cuttv = + A- 

In the long distance limit a 3> /?/ (47r) one has to keep only m = term 
in (EH]) . Thus the free energy of the plates inside a cylinder in the high 
temperature limit is equal to: 

fl)\a>p/(4n) = E ~ ex P(~ 2a W 



kD 



+ ^E ln(l-exp(-2aA iJV )). (36) 



One can check that the limit Axdcl <C 1, Xincl C 1 in ( 1361) immediately 
yields the known high temperature result for two parallel perfectly conducting 
plates separated by a distance a (see eq.( l5Tj) for details). 

If the conditions Ai£>a 3> 1, A^a ^> 1 are satisfied in addition to the con- 
dition a 3> /? / (47r) then the leading asymptotic behaviour of the free energy 
can be expressed via the lowest positive eigenvalues of the two dimensional 
Dirichlet and Neumann problems Aid and Aiat as follows: 

^W)L>/V(4*), w>i, x 1N a>i ~ "^( e " 2aAlD + e" 2aAlAr ), (37) 

so the force between the two plates in a cylinder is exponentially small for 
large enough distances between them at finite temperatures as well. 
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4 Free energy: short distance behaviour 

It is convenient to apply the technique of the heat kernel and zeta function 
to obtain the short distance behaviour of the free energy (1341) . It can be done 
by noting that if the heat kernel expansion 

+00 

A,; fc=0 

exists (d is a dimension of the Riemannian manifold) then one can write the 
expansion 

Te-^ ~ Y *J} d -Q t^c k (39) 

by making use of the analytical structure of the zeta function. 

The proof can be done as follows. Zeta function can be written in two 
different forms: 

1 p+00 

^ = E f5/2) I * tH °p(-^> (40) 

1 r+00 

^ = T,2r(s) I d^'^M-ViXd. (41) 

It is well known that residues at the poles of the zeta function are related to 
the coefficients c k of the heat kernel expansion ( |38i) : 

c k = ^Res s=d _ fc (C(s)r(s/2)), (42) 



which follows from ( j4T)l) and ( l38l) . The expansion ( |39l) now follows from (jH 
and (PI). 

The free energy ( 1341) can be rewritten as follows: 
T{a, (3) = \ ln(l " exp(-2aAg)) + J] ^ ln(l - exp(-2aA feD )) + 



+ E 7? ln ( X - exp(-2aA^)) - £ -L l n (l - exp(-2aA w )) • 

A (3) ^ x m 



*sJV 



- E W 1 - e *p(- 2aA $)) ( 43 ) 



*IN 



10 



Here the eigenvalues A^ satisfy the equation A^p r (x, y, z) = —\^p r (x, y, z) 
and p r (x,y,z) satisfy Dirichlet boundary conditions at the boundary of the 
manifold [0, (3/2} xM, the eigenvalues A^- satisfy the equation A^p s (x, y, z) = 

(3)2 

~Kn Ps( x , y, z) and p s (x, y, z) satisfy Neumann boundary conditions at the 
boundary of the manifold [0, /3/2] x M, M is an arbitrary cross section of the 

(3) 

cylinder. Here we sum over all nonzero eigenvalues X s! ^. 

The eigenvalues A[^ satisfy the one dimensional Laplace equation with 
Neumann boundary conditions at the boundary of the manifold [0, 0/2]. 
They appear due to the condition \ iN ^ in equation (134")) and our decision 
to sum over all nonzero A^ in (j4Hl) . Thus in the last line of fj4*31) the eigenval- 
ues A^ corresponding to the eigenfunctions p n [x,y,z) = cos(27rnz / (3) (for 
which Ajjv = 0) are effectively subtracted. 

Our strategy is the following: one expands the logarithms in the formula 
(|43|) in series 

, exp(— 2an\i) . . 

In 1 - exp(-2oAi) = - V — ^ 44 

L — ' n 

n=l 

applies the expansion flHUj) to each term in f|4*4l) for short distances a and per- 
forms the sum over n thus getting Riemann zeta function at integer positive 
values. 

It is possible to obtain the coefficients of the heat kernel expansion for 
the operators along the following lines. For the manifold [0, /3/2] x M 
and Dirichlet boundary conditions one can write: 




Here 5* is an area and P is a perimeter of the cross section M, a, is the 
interior angle of each sharp corner at the boundary dM and L aa {^j) is the 
curvature of each boundary smooth section described by the curve jj. 

Thus the important for our purpose Seeley coefficients for the operator 
A^ 3 ) with Dirichlet boundary conditions imposed at the boundary of the 
manifold [0, (3/2] x M can be obtained from ( l4"5i) as the coefficients at specific 
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powers of t in the expansion flHHl) : 




42 = d^Ti («) 

45 = ---^- 117) 



One can check by a direct calculation that = cf^, cf^ = — c^, 
c 2iv = C 2D ^ or nianifold [0, /?/2] x M with Neumann boundary conditions. 
The other needed coefficients for manifold M can also be read off from 



ED: c$ = -c$ = -P/(8Vi), for manifold [0,/?/2]: c$ = 
For a <C /3/(47r) and A^a <C 1, Aijva <C 1 one obtains from (1391) and (T43l) 
the leading terms for the free energy: 

•^'(a,/?)|a</3/(47r), A 1I5 a«l, A 1Jv a«l = ~ "^2~ ^3 ~ ^ C 1 ~ 2 X> + 0(1), (49) 



where 

X 



E^-^+E^/w^, (50) 

The force calculated from (j4"9~l) coincides with the zero temperature force Fc 
in [3], (Eq.7). Thus we prove that in the finite temperature case the leading 
short distance terms are the same as in the zero temperature case. 

In the limit Aldo <C 1, Aiatgi <C 1 one immediately obtains from (1361) the 
high temperature result for two parallel perfectly conducting plates separated 
by a distance a. One expands logarithms in series and uses fl39|) and c d = 
Cqn = S/ (47r) in two dimensions (d = 2) to obtain: 



o>/3/(47r), Ai a<l, A 1Jv a<l 

El exp(— 2an\kD) sr^ 1 



2/3^ n 



^ exp(— 2an\iN) 



a^o f— ' 2/3 ^ n 



^111, , Ch(3) ^ 

n=l > \ / 



which is a well known result P2J [13, HE] . 
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